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Abstract: 

<p>Shape  asymmetries  in  nominally  circular  nanomagnets  provide  a  potential  means  for  vortex 
chirality  control.  Flowever,  in  realistic  arrays  their  effects  are  challenging  to  probe  since  asymmetric 
magnetization  reversal  processes  are  often  averaged  to  include  distributions  over  all  angles.  Flere 
we  investigate  how  shape  asymmetry  influences  the  vortex  reversal  in  arrays  of  sub-micron  edge- 
cut  Co  dots.  We  find  that  the  vortices  can  be  manipulated  to  annihilate  at  particular  sites  under 
different  field  orientations  and  cycle  sequences.  The  vortex  annihilation  field  and  degree  of  chirality 
control  depend  sensitively  on  the  angular  position  of  the  applied  field  relative  to  the  flat  edge  of  the 
dots.  For  small  angles,  the  major  loop  annihilation  field  is  significantly  larger  than  that  found  from 
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Abstract 

Shape  asymmetries  in  nominally  cireular  nanomagnets  provide  a  potential  means  for 
vortex  ehirality  control.  However,  in  realistic  arrays  their  effects  are  challenging  to  probe  since 
asymmetric  magnetization  reversal  processes  are  often  averaged  to  include  distributions  over  all 
angles.  Here  we  investigate  how  shape  asymmetry  influences  the  vortex  reversal  in  arrays  of 
sub-micron  edge-cut  Co  dots.  We  find  that  the  vortices  can  be  manipulated  to  annihilate  at 
particular  sites  under  different  field  orientations  and  cycle  sequences.  The  vortex  annihilation 
field  and  degree  of  chirality  control  depend  sensitively  on  the  angular  position  of  the  applied 
field  relative  to  the  flat  edge  of  the  dots.  For  small  angles,  the  major  loop  annihilation  field  is 
significantly  larger  than  that  found  from  the  half  loop  and  the  vortex  chirality  can  be  well 
controlled.  At  intermediate  angles  the  chirality  control  is  lost  and  an  interesting  crossover  in  the 
annihilation  field  is  found:  the  half  loop  actually  extrudes  outside  of  the  major  loop,  exhibiting  a 
larger  vortex  annihilation  field.  At  large  angles  the  annihilation  fields  along  major  and  half  loops 
become  degenerate. 

PACS  number(s):  75.75.+a,  75,60,Jk,  75.60,Ej,  75.70.Kw 
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I.  Introduction 

Control  over  domain  structures  in  magnetic  nanoelements  is  critical  to  the  understanding 
and  applications  of  such  materials/'"'  In  particular,  magnetic  vortices  in  sub-micron  patterns  have 
gained  considerable  interest  in  recent  years  due  to  their  unique  reversal  mechanisms  and 
potential  applications  in  ultrahigh  density  patterned  magnetic  recording  media/''"  The  vortices 
are  characterized  by  an  in-plane  magnetization  with  clockwise  or  counter-clockwise  chirality  and 
a  central  core  with  out-of-plane  magnetization  (up  or  down  polarity).  The  ability  to  control  the 
different  vortex  configurations  within  a  single  nanomagnet  introduces  alternative  data  storage 
possibilities."’  Typically  in  circular  dots  the  vortex  annihilation  field  is  independent  of  where 
the  vortex  core  is  expelled  from  the  dot.  In  realistic  assemblies  of  dots,  variations  in  dot  shape, 
size,  and  intrinsic  anisotropy  inevitably  exist  and  affect  the  reversal  processes.  In  particular, 
dot  asymmetry  has  been  shown  to  lift  the  degeneracy  in  vortex  chirality,  therefore  providing  a 
means  for  chirality  control. Recently,  in  studies  of  arrays  of  67  nm  Fe  nanodots  we  have 
found  distinct  annihilation  fields  depending  on  which  side  of  the  dot  annihilation  occurred.'^’ 
The  primary  cause  is  the  asymmetry  in  the  dot  shape  due  to  deviations  from  perfect  circles.  In 
these  Fe  nanodots,  as  is  typical  in  nanomagnet  arrays,  the  slight  shape  asymmetry  is  randomly 
distributed,  leading  to  asymmetric  reversal  in  individual  dots  but  overall  isotropic  behavior  in  the 
array.  A  key  question  is  how  the  asymmetry  influences  the  vortex  reversal  when  its  orientation  is 
varied  relative  to  the  applied  field.  For  example,  how  does  the  vortex  nucleation/annihilation 
field  change  and  is  the  chirality  control  always  maintained?  In  this  work  we  explore  these  issues 
in  arrays  of  asymmetric  Co  dots  where  the  circular  shape  in  all  the  dots  has  been  broken  in  the 
same  fashion  in  order  to  gain  insight  into  the  reversal  processes. 
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II,  Experiments 

Polycrystalline  arrays  of  Co  dots  were  fabrieated  on  naturally  oxidized  Si  substrates  with 
standard  eleetron  beam  lithography  and  lift  off  teehniques,  in  eonjunetion  with  magnetron 
sputtering.  A  seanning  eleetron  mieroseopy  (SEM)  image  is  shown  in  Fig.  1(a).  The  edge-eut 
dots  are  nominally  40  nm  in  thiekness  and  685  nm  in  diameter.  They  form  a  square  array  over  a 
100x100  pm  area,  with  a  eenter-to-eenter  separation  of  900  nm.  The  pattern  used  to  ereate  eaeh 
dot  during  e-beam  writing  is  a  regular  dodeeagon  with  three  of  the  sides  removed  [Fig.  1(a) 
inset],  thus  ereating  an  asymmetry.  Photo-emission  eleetron  mieroseopy  (PEEM)  studies  were 
earned  out  at  beamline  11.0.1  of  the  Advaneed  Eight  Souree  (AES).  For  eomparison,  arrays  of 
eireular  Co  dots  are  also  examined,  40  nm  in  thiekness  and  870  nm  in  diameter,  with  a  eenter-to 
eenter-spaeing  of  1  pm.  Remanent  state  images  of  both  types  of  dots  are  shown  in  Figs.  1(b)  and 
l(o),  after  saturating  the  dots  to  the  right.  All  the  dots  reverse  the  magnetization  via  a  single 
vortex.  Most  of  the  edge-eut  dots  exhibit  a  eounter-eloekwise  ehirality  [Fig.  1(b)  inset], while 
the  eireular  dots  show  a  random  distribution  of  ehirality. 

Magnetie  hysteresis  loops  were  measured  at  room  temperature  using  the  magneto-optieal 
Kerr  effeet  (MOKE)  on  a  Durham  Magneto  Opties  NanoMOKE2  magnetometer.^  The  beam 
was  foeused  to  about  a  30  pm  diameter  spot  size,  eapturing  the  average  reversal  behavior  of  ~10 
dots.  The  direetion  of  an  in-plane  applied  magnetie  field  was  varied  relative  to  the  fiat  edge  of 
the  dots,  whieh  is  defined  as  0°.  At  eaeh  angle,  major  loops  were  measured  between  ±650  Oe 
and  half  loops  were  measured  over  650  Oe  -  0  Oe  -  650  Oe,  both  with  a  field  spaeing  of  2-4  Oe. 
At  a  field  sweep  rate  of  11  Hz,  typieally  ~10  loops  were  averaged  to  obtain  a  single  hysteresis 
eurve.  Following  prior  proeedures,  the  first-order  reversal  curve  (FORC)  method  was 
employed  for  seleeted  angles.  After  positively  saturating  the  sample  the  applied  field  was 
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reduced  to  a  given  reversal  field  Hr,  the  magnetization  M  was  then  measured  back  to  positive 
saturation  thereby  tracing  out  a  FORC.  This  process  was  repeated  for  more  negative  reversal 
fields  until  negative  saturation  is  reached.  A  mixed  second  order  derivative  of  the  magnetization 

M  {H,  Hr)  was  used  to  generate  the  FORC  distribution  p  =  -d^M{H,Hj^)/2dHdHj^,  which 

21 

capture  the  irreversible  vortex  nucleation/annihilation  events. 

The  experimental  results  were  also  compared  with  micromagnetic  simulations  performed 
with  the  OOMMF  code.  ’  Material  parameters  suitable  for  these  polycrystalline  Co  dots  were 
used  (saturation  magnetization  Ms  =1.4x10^  A/m  and  exchange  stiffness  A  =  3x10'*^  J/m).^"^  A 
weak  residual  uniaxial  anisotropy  of  Ai=l.lxl0"^  J/m^  with  an  easy  axis  along  the  fiat  edge  of  the 
dot  was  found  experimentally  and  included  in  the  simulations.  A  SEM  image  of  an  actual  dot 
was  used  to  construct  the  simulated  dot,  as  shown  in  Fig.  4(c),  thus  reproducing  the  rounded 
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edges  from  the  lithography  process. 

III.  Results  and  Discussions 

The  magnetic  hysteresis  loops  of  the  Co  dots  are  similar  at  different  angles  -  all  have  the 
characteristic  pinched  shape  with  zero  remanence  and  abrupt  magnetization  jumps  that 
correspond  to  the  vortex  nucleation  and  annihilation.^  Representative  sets  of  major  and  half  loops 
are  shown  in  Fig.  2  for  three  angles.  For  clarity,  only  the  first  quadrant  is  shown  [one  full  loop  is 
shown  as  the  outer  boundary  of  FORC’s  in  Fig.  3(a)].  At  0°,  the  vortex  annihilation  along  the 
half  loop,  marked  by  the  abrupt  magnetization  jump  to  positive  saturation,  occurs  much  sooner 
than  that  along  the  major  loop  [Fig.  2(a)].  Surprisingly,  at  60°  [Fig.  2(b)]  the  opposite  behavior 
is  observed  where  the  annihilation  along  the  half  loop  occurs  noticeably  later  than  that  along  the 
major  loop.  This  leads  to  an  unusual  behavior  where  the  half  loop,  a  particular  minor  loop. 
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extrudes  outside  of  the  major  loop.  This  unusual  pattern  is  a  direet  eonsequenee  of  the  dot 
asymmetry  and  the  changing  energy  landscape  during  field  cycling.  The  vortex  annihilation 
along  the  half  loop  faces  a  higher  energy  barrier  than  along  the  full  loop.  Finally,  at  90°  [Fig. 
2(c)]  no  discernable  difference  in  vortex  annihilation  is  observed.  This  trend  was  also 
qualitatively  reproduced  in  simulated  major  and  half  loops,  Fig.  2  (d-f),  for  the  same 
representative  angles.  A  comparison  of  the  calculated  micromagnetic  energies  reveals  that  the 
demagnetization  and  Zeeman  energies  play  a  dominant  role  during  vortex  nucleation  and 
annihilation,  while  the  weak  anisotropy  is  the  least  significant. 

A  related  asymmetry-driven  vortex  annihilation  field  was  previously  inferred  in  arrays  of 
67  nm  Fe  nanodots. Vortex  annihilation  fields  along  first  order  reversal  curves  also  showed  a 
crossover  as  vortices  were  annihilated  from  opposite  sides  of  the  dots.  Consequently,  a 
pronounced  negative-positive-negative  trio  of  features  were  found  at  the  lower  right  corner  of  the 
FORC  distribution  p{H,  where  a  positive  peak  in  p  was  accompanied  by  two  adjacent 

negative  valleys,  highlighting  the  effect  of  shape  asymmetry.  Flowever,  the  exact  angular 
dependence  of  the  annihilation  field  could  not  be  resolved  due  to  the  random  distribution  of  the 
shape  asymmetry  in  the  Fe  nanodots.  For  comparison,  we  have  carried  out  FORC  analysis  on  the 
Co  edge-cut  dots  discussed  here. 

At  0“  the  FORC  distribution  is  characteristic  of  reversal  via  a  vortex  state. As  shown 
in  Fig.  3(b),  a  first  prominent  peak  highlighted  by  region  1  at  (//,  //«)  ~  (430  Oe,  220  Oe) 
corresponds  to  the  annihilation  of  vortices  from  the  flat  side  of  the  dot  and  essentially  maps  out 
the  irreversible  processes  along  the  half  loop,  whose  ascending  branch  is  simply  a  FORC  with 
Hr=0.  a  second  peak  highlighted  by  region  2  at  (-220  Oe,  -450  Oe)  corresponds  to  the 
nucleation  of  vortices  from  negative  saturation.  Highlighted  in  region  3  at  ~  (400  Oe,  -  450  Oe) 
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is  a  pronounced  negative-positive  pair  of  features,  unlike  the  aforementioned  trio  of  features 
observed  in  Fe  nanodots/^’  Note  that  along  sueeessive  FORC’s  with  more  and  more  negative 
reversal  fields,  the  vortex  annihilation  field  approaehing  positive  saturation  moves  progressively 
higher;  over  the  applied  field  range  of  350  Oe  -  500  Oe,  the  slope  of  the  FORC’s  first  deereases 
and  then  inereases,  leading  to  respeetively  negative  and  positive  values  of  p.  The  behavior  has 
also  been  previously  observed  in  simulated  FORC’s  on  edge-eut  Fe  dots  and  is  due  to  the 
differenee  in  annihilation  fields,'^  whieh  for  instanee  ean  be  easily  observed  between  major  and 
half  loops.  This  interesting  annihilation  behavior  beeomes  elear  in  the  FORC  and  half  loop 
analysis  but  is  hidden  when  analyzing  the  major  hysteresis  loops  alone.  At  60°,  a  family  of 
FORC’s  is  shown  in  Fig.  3(a).  The  eorresponding  FORC  distribution  [Fig.  3(o)]  resembles  that  at 
0°,  exeept  for  a  weak  negative-positive-negative  trio  of  features  in  region  3  [Fig.  3(o)  inset].  This 
set  of  features,  whieh  is  almost  identieal  to  that  found  in  the  Fe  dots  diseussed  earlier  and  is 
eaused  by  some  of  the  FORC’s  extruding  outside  of  the  major  loop.  The  weak  intensity  is  a 
manifestation  of  the  small  differenees  in  annihilation  fields  along  sueeessive  FORC’s.  At  90°,  the 
two  main  peaks  remain  in  the  FORC  distribution  [Fig.  3(d)];  however,  the  feature  in  region  3  has 
faded  away  sinee  onee  vortiees  have  nueleated  the  annihilation  field  along  subsequent  FORC’s 
remains  the  same  (the  traee  amount  of  a  residual  feature  is  due  to  the  small  variations  in  the 
array).  In  the  previously  studied  Fe  dots,  where  all  shape  asymmetries  were  randomly 
distributed  and  averaged  over,^^’  the  three  reversal  behaviors  typified  by  the  FORC  diagrams 
in  Figs.  3(b-d)  would  all  eontribute  to  the  observed  negative-positive-negative  trio  of  features. 
These  results  also  demonstrate  that  shape  asymmetry  has  a  distinet  effeet  on  the  vortex 
annihilation  field,  depending  on  the  field  eyele  sequenee,  and  ean  be  turned  on  and  off  by 
varying  the  angular  positions. 
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The  annihilation  field  along  major  and  half  loops  ean  be  determined  quantitatively  from 
the  field  at  whieh  the  magnetization  jumps  abruptly,  i.e.,  where  the  M-H  eurve  has  a  maximum 
slope.  The  angular  dependence  of  the  annihilation  field  extracted  from  the  derivative  of  the 
measured  loops  is  shown  in  Fig.  4(a).  Three  distinct  regions,  represented  by  the  vortex 
annihilation  behavior  shown  in  Fig.  2,  are  found:  at  low  angles  the  annihilation  field  along  the 
half  loop  is  significantly  smaller  than  that  along  the  major  loop;  at  intermediate  angles, 
especially  over  55“-65“,  the  half  loop  annihilation  field  is  slightly  larger;  at  even  higher  angles, 
approaching  90°,  the  two  annihilation  fields  converge.  The  resultant  angular  dependence^^  has  a 
crossover  region,  roughly  corresponding  to  the  angular  positions  when  the  applied  field  passes 
through  the  corners  of  the  fiat  edges  of  the  dots. 

Micromagnetic  simulations  also  show  the  three  distinct  regions  in  the  angular 
dependence,  as  shown  in  Fig.  4(b),  qualitatively  reproducing  the  measured  results.  Since  the 
vortex  nucleation  along  both  major  and  half  loops  is  the  same,  the  difference  in  annihilation  field 
is  better  illustrated  by  examining  the  vortex  core  annihilation  sites  found  from  simulations,  as 
shown  in  Fig.  4(c).  At  small  simulated  angles  (<30°),  the  vortex  annihilation  site  is  well  defined 
due  to  chirality  control  achieved  with  an  asymmetric  dot  [Fig.  1(b)]. During  reversal  from 
saturation  the  magnetization  preferentially  buckles  towards  the  flat  edge  of  the  dot,  assisted  by 
the  demagnetizing  field.  A  vortex  is  nucleated  from  the  fiat  edge  of  the  dot,  and  subsequently 
annihilated  from  the  rounded  edge  of  the  dot  along  a  major  loop.  However,  if  the  field  sweep  is 
stopped  at  zero  field  and  reversed  towards  positive  saturation  (i.e.  tracing  out  a  half  loop)  the 
vortex  core  must  annihilate  from  the  fiat  edge  of  the  dot.  For  angles  larger  than  30°,  simulations 
qualitatively  reproduce  the  crossover  region  where  the  half  loop  annihilation  field  is  larger  than 
the  major  loop  value.  We  find  that  the  chirality  control  is  lost  and  the  annihilation  site  is  either  at 
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a  corner  of  the  flat  edge  or  the  more  rounded  edge  of  the  dot.  For  angles  near  45°,  annihilation 
from  the  eorner  (major  loop)  oeeurs  at  a  smaller  field  than  the  rounded  edge  of  the  dot  (half 
loop),  due  primarily  to  exehange  energy  gains,  leading  to  the  crossover  behavior.  Approaehing 
90°  the  dot  asymmetry  no  longer  plays  a  signifieant  role  as  the  vortex  eore  moves  parallel  to  the 
fiat  edge  and  the  reversal  mimies  that  of  a  symmetrie  dot.  The  annihilation  sites  along  major  and 
half  loops  beeome  degenerate.  Note  that  although  these  simulations  are  illustrative  of  the  typieal 
behavior  in  a  single  dot,  some  differenees  in  the  details  from  experimental  results  are  still 
expected  sinee  the  aetual  arrays  of  dots  do  have  finite  variations  in  their  characteristics.^^ 

It  is  worthwhile  to  examine  the  dipolar  interaetions  in  the  array,  whieh  eould  potentially 

o 

lead  to  a  magnetostatieally  indueed  anisotropy.  Novosad  et  al.  studied  square  array  of  permalloy 
dots  and  found  that  as  the  interdot  distanee  decreased  the  nucleation  and  annihilation  fields 
shifted  towards  zero  field  and  the  magnetostatic  interactions  between  dots  began  to  play  an 
important  role.  They  also  showed  that  the  interaetions  redueed  as  the  dot  aspeet  ratio  (dot 
thickness  L!  radius  K)  deereased.  Guslienko  ealeulated  the  fourfold  anisotropy  eonstant  in 
reetangular  arrays  of  dots  as  a  function  of  normalized  interdot  distanee.’  For  our  dots,  the 
normalized  interdot  spacing  d=dlR=  0.63  and  the  aspeet  ratio  is  T/i?=0.1 1,  where  d  is  the  edge-to- 
edge  dot  spaeing.  Aeeording  to  Refs.  7  and  8,  our  dots  are  approaehing  the  non-interaeting 
regime  and  the  magnetostatieally  indueed  anisotropy  due  to  the  array  layout  is  not  appreeiable. 
In  our  OOMMF  simulations  of  3x3  and  5x5  arrays  of  dots,  the  reversal  behaviors  are  nearly 
identical,  and  are  qualitatively  similar  to  those  found  in  a  single  dot,  as  shown  in  Fig.  5.  Only  a 
very  small  reduetion  in  the  nueleation/annihilation  fields  is  observed  in  the  array  simulations. 
Experimentally,  only  in  an  array  of  symmetrie  circular  Co  dots  with  a  much  smaller  normalized 
interdot  spaeing  of  6=0.21,  do  we  find  evidenee  of  the  magnetostatieally  indueed  anisotropy. 
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Therefore  we  eonelude  that  the  present  asymmetrie  dot  arrays  with  ^=0.63  are  largely  non- 
interaeting. 

IV,  Conclusions 

In  summary,  we  have  found  that  in  asymmetric  Co  dots  the  vortex  annihilation  field  and 
degree  of  chirality  control  depend  sensitively  on  the  angular  position  of  the  applied  field  relative 
to  the  flat  edge  of  the  dots.  For  small  angles,  the  vortex  is  more  easily  expelled  from  the  fiat 
edge  of  the  dots  along  a  half  loop  than  from  the  rounded  edge  of  the  dots  along  a  major  loop. 
The  large  difference  between  annihilation  fields  can  be  used  to  identify  the  vortex  chirality.  At 
intermediate  angles  the  chirality  control  is  lost  and  the  opposite  trend  is  observed.  Along  the  half 
loop  vortex  annihilation  from  the  rounded  edge  of  the  dots  is  harder  than  that  from  the  dot 
comers  along  the  major  loop.  Finally,  at  large  angles  approaching  90“  the  dot  asymmetry  is 
effectively  removed  as  the  vortex  core  moves  parallel  to  the  fiat  edge.  Our  results  demonstrate 
an  intrinsic  effect  of  the  shape  asymmetry  and  illustrate  how  the  vortices  can  be  manipulated  to 
annihilate  at  particular  sites  under  certain  field  orientations  and  cycle  sequences. 
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Figure  Captions 


Fig.  1 .  (Color  online)  (a)  SEM  image  of  arrays  of  Co  dots  with  a  horizontal  flat  edge.  The  inset 
shows  the  pattern  used  by  the  e-beam  writer  to  ereate  eaeh  dot  and  the  0°  orientation  of 
the  sample.  PEEM  images  of  a  typieal  portion  of  (b)  the  edge-eut  dots  and  (e)  referenee 
eireular  dots  at  zero  field  after  saturating  the  dots  to  the  right.  All  dots  are  in  the  single¬ 
vortex  state.  The  ehirality  is  eontrolled  in  (b),  as  shown  in  the  inset,  and  random  in  (e). 

Eig.  2.  (Color  online)  (Eeft)  Measured  major  and  half  loops  using  MOKE,  and  (Right)  simulated 
major  and  half  loops  with  the  applied  field  at  (a,  d)  0°,  (b,  e)  60°  and  (c,  f)  90°  relative  to 
the  flat  edge  of  the  dots. 

Eig.  3.  (Color  online)  (a)  Measured  EORC’s  and  (b-d)  EORC  distributions  with  the  applied  field 
at  (b)  0°,  (a,e)  60°  and  (d)  90°  relative  to  the  flat  edge  of  the  dots.  Cireles  in  (b)  highlight 
the  3  regions  of  EORC  features.  Inset  in  (c)  shows  a  zoom-in  view  of  the  negative- 
positive-negative  set  of  features. 

Eig.  4.  (Color  online)  Angular  dependenee  of  the  annihilation  fields  from  both  major  (solid 
squares)  and  half  loops  (open  eircles)  extraeted  from  (a)  MOKE  measurements  and  (b) 
mieromagnetic  simulations.  Error  bars  are  ineluded  to  indieate  the  field  spaeing  of  the 
measurements  and  the  lines  are  guides  to  the  eye.  The  loeations  of  the  simulated  vortex 
eore  annihilation  sites  are  shown  in  (e).  Eines  are  used  to  highlight  the  annihilation  sites 
for  eaeh  angular  position. 

Eig.  5.  (Color  Online)  Simulated  hysteresis  loops  for  a  5x5  array  (blue  triangles),  3x3  array  (red 
open  eircles),  and  single  dot  (black  squares). 
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Fig.  1.  Dumas  et  al. 
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Fig.  2.  Dumas  et  al. 
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Fig.  3.  Dumas  et  al. 
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Fig.  4.  Dumas  et  al. 
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Fig.  5.  Dumas  et  al. 
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